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I. INTRODUCTION 



Motivated mainly by the Bartnik-McKinnon discovery pLj] of globally regular, static, 
spherically symmetric solutions of the Einstein- Yang-Mills (EYM) equations and similar 
results obtained in the EYM and YM theory coupled to a dilaton field 0, i, i, S, H, 
the Euclidean version of the EYM and YM theory coupled to a dilaton were investigated 
recently [1], [9|. While in the four dimensional space-time with Lorentzian signature electric 
part of non-abelian SU (2) gauge field should necessarily vanish for asymptotically flat solu- 
tions 0,1 lllj l situation is changed in Euclidean sector: electric field here plays role similar 
to the Higgs field of the Lorentzian sector 12j and there are nontrivial dyonic type solutions 

a, 0. 

Dilaton field in these investigations was considered to be massless. On the other hand 
there are very strong experimental bounds on dilatonic coupling constant for the massless 
dilaton case and it is thought that at some stage dilaton should obtain potential possibly 
with stable or metastable minimum. Since close to the minimum potential can be well 
approximated by the dilaton mass term it is of certain interest to investigate Euclidean 
solutions in the YM theory coupled to a massive dilaton, which is exactly the aim of 
present Letter. 

The rest of the Letter is organized as follows: in next Section we formulate our model, 
derive equations of motion and asymptotic behaviour of finite action solutions. In Sect. 
Ill we discuss some special solutions of our system. In Sect. IV we present our numerical 
results. Sect. V contains concluding remarks. 



II. THE MODEL 

Starting point of our investigation is the Euclidean version of Yang-Mills theory coupled 
to a dilaton field 0, which is defined by the Euclidean action 

Se = J i^-^d,<pd^^ + J^F.^F^" + Vi<P)j d'x , (1) 

where F^,^ is the non-Abelian gauge field strength, V{(f)) is dilaton field potential and 7 and 
g denote the dilatonic and gauge coupling constants, respectively. 

The simplest choice of a dilaton potential would be to have just a mass term: 

Vi<P) = ^m'<p' , (2) 

where m is dilaton mass. 



A. The Ansatz and equations 



In this study we will restrict ourselves with the SU{2) gauge g roup and will be inter- 



ested in spherically symmetric solutions. Following Witten |l3l . Il4l | the general spherically 



symmetric SU{2) Yang-Mills field can be parameterized as follows: 



W{r) 



Xj J OC Q 



A, 



+ 



(3) 
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with r = \fxl being the radial coordinate, z, j = 1, 2, 3 are spatial indices and a,b = 1, 2, 3 
group indices. 

The reduced action Sred is defined as follows: 



Tmax 



Se = drSred , (4) 







with T = X4 being Euclidean time. Without loss of generality we can put Ai = A2 = 0. 
Then the reduced action reads: 



Sred — I drCymd ; (5) 

with 



where 



1 fl _ U/2\2 1 

^ym = ^(^'^ + + lr\- + W\^) , (7) 

and the prime denotes the derivative with respect to r. 

The equations of motion which follow from the reduced action ([5]) read: 

W" = -2;4.'W'-^1^^ + Wu\ (8) 

u' 2 

u" = -2 270'm' + -^W\ , (9) 

2 , 27e2^<^, „ il-W^f r^u'^ 9 9^ 9 



The rescaling 



<P 1 9 9 

— , r — s> — r , M — s> — M , m —m , (11) 

19 7 7 



removes the gauge and dilatonic coupling constants g and 7, respectively, from the equations 
of motion. We thus set without loosing generality 7 = = 1 in the following. 

B. Asymptotic behaviour 

Close to r = we find a 3-parameter family of solutions regular at the origin, which can 
be parameterized in terms of b, Ui and (po: 



W{r) = 1 - for^ + O(r^) , (12) 
u{r) = Mir- yr2 + 0(r3) , (13) 

0(r) = 00 + (^mVo + 2e2'^" (b^ + ^) y + O^) . (14) 
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The behaviour of the solutions at infinity is more involved. In leading order we have: 

W{r) = Ce-^V'3'=(l + 0(-)) , (15) 

r 

u{r) = f/-^ + 0(l), (16) 

So, at infinity solutions are characterised by three free parameters U^Qe and C. 

Note that in the massless case the shift of dilaton field 0^0 + for any finite value 
can be compensated by a rescaling of r and u as follows: 

r ^re^ , u ue~^ . (18) 

This symmetry is obviously broken by dilaton mass term and for finite action solutions we 
have to choose normalisation 0oo = 4>{r = oo) = 0. 



III. SPECIAL SOLUTIONS 

There are few special solutions of our system, which play an important role. 
The simplest is the vacuum solution with 

W{r) = ±1 , u{r) = , 0(r) = . (19) 

Note here that W{0) = 1 and W{0) = —1 are gauge equivalent and in the Eq. (|T2i) we have 
chosen positive sign as in most of the investigations on the subject. 

In massless case the equations of motion have elegant solution [g], 0, T2| , for which the 
gauge field amplitudes are simply 

W{r) = , u{r) = U , (20) 
with U an arbitrary constant and dilaton field = 0^°^ has logarithmic behaviour 



do) I 



1, 



^ab (r) = -ln(l + -) . (21) 

In massive case this abelian solution is modified as follows: gauge amplitudes W{r) and u{r) 
are the same as in Eq. (120!) and dilaton field = 0^™'' now satisfy massive equation [15] 



2 e2* 



, (22) 



which unfortunately cannot be integrated analytically. Solution of this equation 0^™-* plays 
an important role, while as we discuss below in some limit solution of full system will be 
approaching this abelian solution. 

In the regime, when dilaton decouples, = 0, one finds the BPS monopole solution. Note 
that in fiat space-time there are no radial excitations of BPS monopole [l6| in contrast to 



case when gravity is taken into account [17|, [18 
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IV. NUMERICAL RESULTS 



We integrated the equations of motions fl8l)- f|T0l) using the differential equations solver 



COLSYS [l9j. To construct solutions we used the following boundary conditions: 

W{Q) = 1 , n(0) = , 0(r)'|,=o = (23) 

at the origin and 

W{oo) = , u{oo) = U , 4>{oo) = (24) 

at infinity. 

From numerically found solutions with fixed U and n, the number of nodes of gauge 
field amplitude W, one can extract the values of b, ui and other parameters which enter 
in the asymptotic behaviour Eqs. (|T2|) - f|T^ and Eqs. f|T5|) - f|T7|) . Examples of solutions with 
zero, one and two nodes of the gauge field amplitude W{r) are shown on Fig. [H where we 
plot profiles of W{r),u{r) and 0(r) for U = 1.0 and = 10.0. Typically dilaton field 
grows monotonically from some (negative) value 0o at r = to 0oo = at r = cxd. Electric 
component of the YM field, u(r), also grows monotonically from u{0) = to some u{oo) = U 
and gauge amplitude W{r) starts from its vacuum value W{0) = 1 and for the solution with 
n nodes oscillates n times and then tends to zero according to Eq. f|T5l) . 

In order to interpret our numerical results, we find it convenient to analyze the Yang-Mills 
part of the action density, i.e. the piece Cym, Eq- ©5 appearing in the full action density 
Eq. ([6]). There are two factors which can make Lym large in the neighbourhood of the core 
of the soliton: (i) the case U ^ 1 , (ii) when many nodes appear, the function W and/or 
the ratio (1 — W'^)/r can get large and lead to important contribution to the action density. 

The consequence of these features on the pattern of solution is observed already in the 
massless case [9[. For fixed n and increasing the parameter U, we observe a decreasing of the 
value 00- The factor exp{2(f)) then get very small and decreases considerably the contribution 
to the action of Lym in the region of the core of the soliton (i.e. for r ~ 0). 

In the case of a massive dilaton, the situation is more complex because, increasing too 
much the value of the parameter m also leads to a large term in the action density and there 
appears to be a competition between the Yang- Mills part and the dilaton part. 

Note that in problem under consideration we have three different scales: YMD scale, 
^ymd = ^, 'monopole' scale, Imon = jj and scale associated with the dilaton mass, IdUaton = 
With the rescaling (ITTl) we effectively put lymd = 1, so we still have two scales to vary. It 
turns out that the large behaviour of nodeless solutions and solutions with the nodes are 
different. So we consider these cases separately. 



A. Zero-node solutions 

The result of the competition between the Yang-Mills field and the massive dilaton is 
illustrated by Fig. [21 where the value of 0o is plotted as function of for several values of 
U for lowest branch of solutions with no nodes. 

For small values of U, we see that increasing the mass of the dilaton has the effect to 
increase the value 0o monotonically. As a consequence, for large enough dilaton masses, the 
function (j){r) becomes uniformly small and the corresponding Yang-Mills field converges to 
the BPS monopole, suitably rescaled, since U 1. For U > 1 the scenario is different, 
as suggested by Fig. [2j Indeed, increasing m, we observe that the value of 0o becomes 
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more negative in the region of the origin, so that the large values of Lym are lowered by the 
dilaton factor. While keeping m increasing, we observe that the value 0o reaches a minimum 
for some m = rricr and then increases to reach zero. We would like to point out that our 
numerical analysis becomes rather involved for > 1000. The occurrence of a minimal 
value for 0o can be related to the balance between the Yang-Mills and dilaton pieces of the 
action density. 

The behaviour of the solutions in the large U limit and with m > fixed is illustrated by 
Fig. [3] in the case n = and = 10. Here we reported the Yang-Mills action density and 
the dilaton function (f){r) for three values of U, namely U = 2, 10, 100. We see in particular 
that the function Lym reaches its maximum closer to the origin while increasing U defining 
the radius of the soliton core. Outside this core the dilaton field approaches the profile 0^^'* 
(which is m dependant) in a larger and larger region of space. Only inside the core does the 
dilaton deviates from 0^7'*- Interestingly, our numerical analysis strongly suggest that for 
large U, the relation 

(f)Q = F - \n{U) , with F = const for U ^ oo (25) 

holds, for both cases m = and m ^ 0. In the massless case this relation leads to the 
existence of critical value tier when we transform to solutions to an alternative gauge, say 
0, If where the dilation shift symmetry Eq. 0181) is fixed by demanding 0(0) = 0. The critical 
value Ucr is then determined by exp(F) where F is the constant above. Our result strongly 
suggest that the constant F is independent on m; at least it was checked with a good level 
of accuracy for < < 100. 

This remarkable numerical result shows that the critical value are independent on the 
mass and simplifies considerably the analysis of the domain of the solutions. It is illustrated 
on Fig. m where the value 0o is plotted as a function of the parameter 5 = 0o + ln(f/) for 

= 0, 10, 100. We have for this case F ^ —1.08. 

The values of Euclidean action A = Sred and parameter Qe are plotted as function of 
ln([/) for = 0, 10, 100 for zero nodes solutions on Fig. ([5]). 



B. Solutions with the nodes 

In w? > case the node solutions constructed in [9i] get smoothly deformed by the dilaton 
mass parameter. Keeping U fixed and increasing the mass, we observe that the value of the 
dilaton function at the origin becomes more and more negative. This is due to the fact that, 
the position of the node, say Tq, of the function W, is very close to the origin, so that the 
terms W and (1 — W'^Y/r'^ are quite large in this region and favorize negative values of 0(r) 
for r e [0,ro]. This is illustrated by Fig. [H] where the values 0o are reported as function of 

for U = 0.1 and U = 1.0 (for larger U the pattern is the same). Note that this behaviour 
of 00 of Fig. [6] is different from the zero- node case Fig. [2] and this difference is explained by 
absence of radial excitations of BPS monopole in flat space-time (l6| . 

The large [/-limit of one node solutions was also investigated for several values of m and 
is illustrated on Fig. [71 We see in particular on this figure that the property (125|) also hold 
for one-node solution. Here we find F ^ —4.70. 

The values of Euclidean action A = Sred and parameter are plotted as function of 
ln{U) for = 0, 10, 100 for one nodes solutions on Fig. ([H]). 
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V. CONCLUDING REMARKS 



In the present study we investigated 0(3)— symmetric solutions of four dimensional Eu- 
clidean YMD theory and found branches of dyonic type solutions for any values of dilaton 
mass m. They can also be viewed as static (vortex type) solutions of 4 + 1 dimensional 
theory with Lorentzian signature. Note that in pure Euclidean YM theory solutions with 
higher, 0(4) symmetry are well known instantons. Simple scaling arguments show that there 
should not be any finite action 0(4) symmetric solutions in combined YM-dilaton theory. 

The electric component of YM field in Euclidean theory plays role similar to the Higgs 
field and correspondingly present model is very similar to the static, spherically symmetric 
sector of the YM-Higgs model investigated by Forgacs and Gyiiriisi (FG) [l5|. The main 
difference with the FG model is in the way the dilaton field couples to the the Higgs part of 
the lagrangian. 

It was shown [23] that static spherically symmetric solutions of YM dilaton theory found 
in [s], [g], 0] have sphaleronic nature similar to their EYM "relatives" 21, 22, 2^, namely 



they have unstable mode(s), half-integer topological charge and there are fermionic zero 
modes in the background of these solutions. Interpretation of Euclidean solutions discussed 
in the present Letter depend in particular on number of negative modes in spectrum of small 
perturbations about them. We plan to address this topic in further work. 

Another interestingtopic for further study is the investigation of Euclidean generalisations 
of axially symmetric 24 1 and platonic |2^ solutions of YMD theory. 
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FIG. 1: Examples of solutions with zero, one and two nodes of the gauge field amplitude W{r) 
plotted for parameter values U = 1 and = 10. 
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FIG. 2: Parameter (/iq is plotted as a function of for different values of U for solutions with zero 
nodes. 
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FIG. 3: The Yang-Mills energy density and the dilaton function are plotted as functions of r for 
= 10.0 and U = 2.0, 10, 100 respectively by the solid, dashed and dotted lines for zero nodes 
solutions. 
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FIG. 4: The value is plotted as function of ln(?7) and of the quantity 5 = (/>o + ln(C/) for 
= 0, 10, 100 for zero nodes solutions. 
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FIG. 5: The value of Euclidean action A = Sred and parameter Qe are plotted as function of \n{U) 
for = 0, 10, 100 for zero nodes solutions. 
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FIG. 6: The value of (j)o is plotted as function of for one node solutions for U = 0.1, 1.0. 
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FIG. 7: The value (pQ is plotted as function of the quantity 6 = (j)o + ln(U) for one node solutions 
for = 0, 10, 100. 
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FIG. 8: The value of Euclidean action A = Sred and parameter Qe are plotted as function of \n{U) 
for = 0, 10, 100 for one nodes solutions. 
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